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Abstract 

A class of spectral problems with a hidden Lie-algebraic structure is considered. 
We define a duality transformation which maps the spectrum of one quasi-exactly 
solvable (QES) periodic potential to that of another QES periodic potential. The 
self-dual point of this transformation corresponds to the energy-reflection symmetry 
found previously for certain QES systems. The duality transformation interchanges 
bands at the bottom (top) of the spectrum of one potential with gaps at the top 
(bottom) of the spectrum of the other, dual, potential. Thus, the duality transfor- 
mation provides an exact mapping between the weak coupling (perturbative) and 
semiclassical (nonperturbative) sectors. 



1 Introduction 



Quasi-exactly solvable (QES) systems are those for which some finite portion of 
the energy spectrum can be found exactly using algebraic means |l|, H [3|, f|, ||. A 
positive integer parameter J characterizes the 'size' of this exact portion of the spec- 
trum. This integer J has both an algebraic significance, related to the dimension 
of a representation, and a geometrical significance, the genus of a Riemann surface 
associated with the spectrum. In ||, the large J limit was identified as a semiclassi- 
cal limit useful for studying the top of the quasi-exact spectrum. It was found that 
remarkable factorizations reduce the semiclassical calculation to simple integrals, 
leading to a straightforward asymptotic series representation for the highest QES 
energy eigenvalue. The notion of energy-reflection (ER) symmetry was introduced 
and analyzed in @: for certain QES systems the QES portion of the spectrum is 
symmetric under the energy reflection E — > —E. This means that for a system 
with ER symmetry, there is a precise connection between the top of the QES spec- 
trum and the bottom of the spectrum. Coupled with the semiclassical large J limit, 
the ER symmetry relates semiclassical (nonperturbative) methods with perturbative 
methods 0. In this present paper, we show that for a class of periodic QES poten- 
tials the ER symmetry is in fact the fixed point (self-dual point) of a more general 
duality transformation. The duality transformation we consider is analogous to, but 
importantly different from, a similar transformation suggested in Ref. || relating 
different QES potentials. 

For periodic potentials, usually, the Lie-algebraic QES sector consists of band 
boundaries which correspond either to periodic or antiperiodic boundary condi- 
tions [J Thus, only one boundary of each band is algebraic. However, for the elliptic 
periodic potentials to be considered below, both boundaries are algebraically cal- 
culable. For given J, there are 2 J + 1 energy eigenvalues in the algebraic sector; 
there are J + 1 allowed bands, with the highest band being open (it stretches up 
to E — oo). The 2 J + 1 algebraic energy levels give all band edges, so that in 
the problem at hand all band edges are algebraically calculable. This fact is well- 
known in the literature. In this paper the main emphasis is on physical aspects 
which have not been investigated so far. When the parameter J becomes large, the 
theory becomes weakly coupled at the bottom of the spectrum and semiclassical at 
the top. We develop this idea — using 1/J as an expansion parameter — in the 
application to the elliptic periodic potentials. We show that the duality between 
weak coupling and semiclassical expansions applies not just to the asymptotic series 
for the locations of the bands and gaps, but also to the exponentially small widths 
of bands and gaps. We observe that the 2J+ 1-level algebraic sector splits into four 
completely disconnected algebraic subsectors. Therefore, instead of diagonalizing a 
(2J + 1) x (2 J + 1) matrix, it is sufficient to diagonalize four matrices which are 
roughly four times smaller in their linear sizes. In addition to the eigenvalues, we 

lr The simplest periodic QES problems were discussed in the literature long ago in particular, 
the Lame system was the subject of investigation in Refs. 10, [LT|. 
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determine, purely algebraically, the eigenf unctions corresponding to the band edges. 
For even J there are J + 1 periodic eigenfunctions and J antiperiodic ones. For odd 
J there are J periodic eigenfunctions and J + 1 antiperiodic ones. 

In Sect. ^]we define our duality and self-duality transformation, and review some 
basic facts about the periodic Lame equation. The duality transformation has both 
an algebraic and geometric interpretation. We show how the duality transformation 
extends the algebraic ER symmetry formalism introduced in . Section |3] is devoted 
to a detailed consideration of four distinct s/(2)-based algebraizations relevant to the 
algebraic determination of the band edges and the corresponding wave functions. We 
explain why the cases of even and odd J must be treated separately and how periodic 
and antiperiodic wave functions emerge. The periodicity versus antiperiodicity is 
due to properties of the quasiphases in these two cases. In Sect. [| we study several 
explicit examples in detail. Section || is devoted to an analysis of the self-dual case. 
In Sect. |5] we describe the duality relation between perturbative and nonperturbative 
approximation techniques, as applied to the computation of both the locations and 
widths of bands and gaps. Finally, Sect. |j] contains a summary and some comments 
regarding further generalizations. 



2 Duality and self-duality in Lame model 

In this section we introduce the notions of duality and self-duality for the following 
one-dimensional quasi-exactly solvable (QES) Lame equation: 

j-^ + J(J + 1) z,sn 2 (<^) - l -J{J + 1) | =E^) . (1) 



Here sn(0|z/) is the doubly-periodic Jacobi elliptic function [|T^, [L3|, the coordinate 



G R 1 , and E denotes the energy eigenvalue. The real elliptic parameter v lies in 
the range < v < 1. The potential in (|l|) has period 2K(u), where 

K{u) = / r 2 

Jo Vl-vsm 2 9 

is the elliptic quarter period. Note that the parameter v controls the period of the 
potential, as well as its strength. As v — > 1, the period 1K[y) diverges logarithmi- 
cally, 2K{y) ~ ln(jj^j), while as v — >• 0, the period tends to a nonzero constant: 
IKiv) — > it. In the Lame equation ([]]), the parameter J is a positive integer (for 
non- integer J, the problem is not QES). This parameter J controls the depth of 
the wells of the potential; the significance of the constant subtraction — ^ J(J + 1) 
will become clear below. As an illustration, we have plotted in Fig. [TJ the potential 
energy in the Schrodinger equation (|l|) for two values of z/, namely, v = 0.95 and 
v = 0.05. 

It is a classic result that the Lame equation ([[J) has bounded solutions ^(0) 
with an energy spectrum consisting of exactly J bands, plus a continuum band 
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Figure 1: The potential energy in the Schrodinger equation (]T|) as a function of 4>. 
The solid curve corresponds to elliptic parameter v = 0.95, for which the period is 
2K(0.95) ~ 5.82. The dashed curve corresponds to v = 0.05, for which the period 
is 2iT(0.05) ~ 3.18. In each case, J = 5. Note how different the two potentials are; 
and yet, their spectra are related by the duality transformation (^). 



[T2]| . It is the simplest example of a "finite-gap" JTj, [15| model |, there being 
just a finite number, J, of "gaps", or "excluded bands" in the spectrum. This 
should be contrasted with the fact that a generic periodic potential has an infinite 
sequence of gaps in its spectrum [M. We label the band edge energies by Ei, with 
I — 1, 2, ... , (2 J + 1). Thus, the energy regions, E 2 i-i < E < E 2 i, and E > E 2 i+i, 
are the allowed bands ("bands" for short), while the regions, E 2 i < E < E 2 i+i, and 
E < E\, are the exclusion bands ("gaps" for short). The wave functions ^(<f>) at the 
band edges are either periodic or antiperiodic functions, with period 2K(i/), while 
within the bands the wave functions are quasi-periodic Bloch functions, which can 
be expressed in terms of theta functions 



12]. 



Another important classic result [IT], |Tj| [TJ| concerning the Lame model ([I]) is 
that the band edge energies Ei, for I = 1, . . . , 2 J + 1, are simply the eigenvalues of 
the (2 J + 1) x (2 J + 1) matrix 



H = Ji + vJi 



DI 



(3) 



where J x and J y are su(2) generators in a spin J representation and I is the unit 

2 Much of our discussion generalizes to other finite gap potentials, but we concentrate here on 
the Lame system for the sake of pedagogical definiteness. 
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matrix. Thus the Lame band edge spectrum is algebraic, requiring only the finding 
of the eigenvalues of the finite dimensional matrix H in ([3]). When v = or 1, it is 
straightforward to write down the eigenvalues of H . But for < v < 1, the matrix 
J 2 + vJy cannot be simply related to the su{2) Casimir, J 2 , and so the eigenvalues 
of H are, in fact, nontrivial functions of v. For example, for J = 1 and J = 2, the 
eigenvalues of H are: 



E 1 


= -l + i/ , 




E2 


= 0, 




E 3 


= v\ 




Ei 


= — l + 2i/ — 2Vl — 


1/ + z/ 2 


Ei 


= -2 + i/, 




E 3 


= -2 + 4z/, 




£4 


= 1 + z/, 




£5 


= -1 + 2Z/ + 2V1- 


Z/ + Z/ 2 



The band (gap) structure in the second case, J = 2, is presented in Fig. 0. 

ENERGY BANDS, J=2 

3 r 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
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Figure 2: The energy bands for the Lame system (fj), as a function of v. The plot 
given is for J = 2. The first (the lowest) allowed band (unshaded) is of the p-ap 
type, while the second one of the ap-p type, see text. 

The shaded areas on this plot (and other similar plots in Figs. [|, [| and ||) are 
the forbidden bands (gaps), while the unshaded areas are the allowed bands. The 
allowed bands can be of two types: p-ap and ap-p. In the allowed bands of the first 



4 



type the lower boundary of the band is determined by a periodic wave function, 
while the upper boundary by an anti-periodic one. In the allowed bands of the ap-p 
type, the band edge structure is reversed: the lower boundary corresponds to an 
anti-periodic wave function. The last allowed band has no upper boundary — it 
stretches up to infinitely high energies. 

In ||20|| , the algebraic form (|3]) of the Lame system was exploited to provide a 
precise analytic test of the instanton approximation, as is discussed further in Sect. ||. 

In this paper we study a special duality of the spectrum of the Lame system ([[]), 
which can be stated succinctly as: 

E[v\ = -E[l-v\ (5) 

where E[v] denotes the spectrum for the potential with elliptic parameter v. That 
is, the spectrum of the Lame system (]!]), with elliptic parameter u, is the energy 
reflection of the spectrum of the Lame system with the dual elliptic parameter 1 — u. 
In particular, for the band edge energies, Ei, which are the eigenvalues of the finite 
dimensional matrix H in (|3|), this means that (for Z = 1,2,...,2J+1) 

Ei[v\ = -Eij+2-i [1 - u] (6) 

This duality can be seen directly in the eigenvalues of the J = 1 and J = 2 examples 
in (f|) . The proof for the band edge energies is a trivial consequence of the algebraic 
realization (^), since 

H[u) = J 2 x + vJ 2 - l -J{J+l)\ 
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(J2 + (1 _^J2)_ Ij(J+l)I 



(7) 



Noting that \J 2 Z + (1 — v) J 2 ] has the same eigenvalues as [J 2 + (1 — v)Jy], the duality 
result (|^) follows. It is instructive to see this duality in graphical form, by looking 
at Figs. |2|, [|, [5] and ^|, which show the band spectra, as a function of u, for various 
different values of J. In each case, the transformation v — *■ 1 — u, together with 
the energy reflection E — ► —E, interchanges the shaded regions (the gaps) with the 
unshaded regions (the bands). 

The fixed point, v — |, is the "self-dual" point, where the system maps onto 
itself, and the energy spectrum has an energy reflection (ER) symmetry, as was 
studied in 0. In Sects. 3 and 4 we discuss in detail how this Lame model fits into 
the energy reflection symmetry classification of Shifman and Turbiner. 

In fact, the duality relation (|5]) applies to the entire spectrum, not just the band 
edges (^|). The easiest way to see this is to realize that the duality relation is a simple 
consequence of Jacobi's imaginary transformation, applied to the Lame equation ([[]). 
The original Lame system (P was defined for real coordinate 0, but we can extend 
into the complex plane. The potential is doubly periodic, with real period 2K 
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Figure 3: The period parallelogram for sn 2 (0|i/) is the rectangle with the sides 
2K(i>) and 2K'[y) = 2K(1 — v). The original Lame equation ([!]) is defined for real 
<ft, on the horizontal thick line. The dual Lame equation is defined for <fi on the 
vertical thick line. The parametric coordinate rotation is (||), and leads to the dual 
Lame equation ([T0|). 



and imaginary period 2K' = 2K(1 — v), as is shown in Fig. |3|. Consider rotating 
the coordinate off the real axis by the transformation 

(jj — — K — i K') , (8) 

so that real can be considered as a parametric coordinate for <p running along a 
vertical line through the point K (the vertical thick line in Fig. |3p. Note that along 
this vertical line the potential sn 2 (</>|i/) is again real. Moreover, using the properties 
of the Jacobi elliptic functions, |12|, [L3| , we get 

usn 2 (K + iK' -z0» = 1 - (1 - z/)sn 2 (0'|l - v) . (9) 

The Lame equation on the real axis maps to another, dual, Lame equation on the 
vertical thick line (see Fig. [3D. In other words, the change of variables @ transforms 
the Lame equation (0) into: 

id 2 1 1 

+ J(J + l)(l-z/)sn 2 (0 / |l-z/)--J(J+l)U(0 / ) = - J B^(0 / )> (10) 



/2 
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So solutions of the Lame equation ([!]) are mapped to solutions of the dual equation 
([10|), v — > 1 — u, and with a sign reflected energy eigenvalue: E — > —E. Returning to 
the complex plane diagram in Fig. |3], the transformation (|S]) clearly interchanges 
the real and imaginary periods. In the self-dual case, u — ^, these two periods are 
equal, so the periodic lattice shown in Fig. |3| is a square lattice. This is the so-called 
lemniscate case []. 

It is important to distinguish our duality transformation from another duality 
transformation considered in || , for which the duality transformation was <p' — i (p, 
which maps the real axis onto the imaginary axis. The Lame potential is also real 
along the imaginary axis. However, along the imaginary axis there is a singularity 
at (j) = i K' . The transformation = icf) maps the original Lame potential to a very 
different potential; so there is no notion of "self-duality" for this transformation. 

To see why bands and gaps are interchanged under our duality transformation 
we recall that the two independent solutions of the original Lame equation (|T]) 
can be written as products of theta functions 



./ 



*±($ = n 



± a.j) 



exp {^(j)Z{otj)) 



[ID 



where the parameters ctj, for j = 1, . . . , J, satisfy a complicated set of J nonlinear 
constraints together with the energy E, and Z(a) is the Jacobi zeta function [O 



Under the change of variables (H) these theta functions map into the same theta 
functions, but with dual elliptic parameter. However, they map from bounded to 
unbounded solutions (and vice versa), because of the "i" factor appearing in (§). 
Thus, the bands and gaps become interchanged. This interpretation of the duality 
transformation in terms of the Jacobi imaginary transformation will be important 
in Sect. | when we discuss WKB techniques. 

Away from the self-dual point, v = |, the duality transformation relates the 
spectrum of one elliptic potential with the spectrum of a different elliptic potential. 
For example, the two potentials plotted in Fig. [I], for v = 0.95 and v = 0.05, 
have energy spectra that are inversions E — > —E of one another. It is striking 
how different these two potentials are, even though their spectra are related by this 
simple energy-reflection symmetry. In the semiclassical large J limit, the duality 
symmetry (^) connects properties of low-lying bands of one potential with high- 
lying gaps of the dual potential. At large J, the barriers between neighboring 
potential wells become high, so that tunneling effects are suppressed. Thus, low- 
lying bands are exponentially narrow, so that we can study their location and their 
width. Duality relates these results to the location and width of high-lying gaps for 
the dual potential. 

There is another interesting semiclassical limit that can be studied for the Lame 
models. Since the elliptic parameter v controls the period, 2K(u), of the potential, 
the duality transformation relates a system with v — > 1, for which the individual 

3 Note that K{l/2) = t where r is the parameter defined in Eq. 
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wells are far separated (recall that K{y) — > oo as v — > 1), with a dual system having 
elliptic parameter v 1 — 1 — v — > 0, for which the wells are shallow and close together. 
In the first case tunneling is suppressed since the neighboring wells are far apart, 
so that semiclassical techniques are appropriate. On the other hand, in the dual 
system the potential becomes weak so that perturbative techniques can be used. 
The duality relation provides a direct mapping between these different approximate 
methods, as is studied in Section ||. 



3 Algebraic approach to QES spectral problem 
with elliptic potentials — generalities 

In this section we show how the duality and self-duality properties are character- 
ized using the algebraic language of QES systems. In accordance with the general 
strategy of generating sZ(2)-based QES spectral problems [0], we introduce a new 
variable 7/(0) satisfying the condition 

(jt) =4(1-1017 [(l-iO + i/irf. (12) 

The solution of this condition appropriate for our purposes is 

rj = l-sn 2 (0|z/). (13) 

For real <fi the variable r/ varies between and 1. As we will see shortly, of relevance 
are functions ^{rj) which are constructed of polynomials of r\ and may have square 
root singularities of the type ^ff\ or a/1 — rj; these functions are either periodic or 
antiperiodic in 0, and can be used to construct band-edge wave functions. 

In terms of the variable i], three s/(2) generators have a differential realization 

T+ = 2jr ] -r ] 2 d v , T° = -j rj + V d v , T~ = d n (14) 

where 

[T + ,T"] = 2T°, [T + ,T°] = -T + , [T~,T°}=T-, (15) 



and = j^. If j is semi- integer, the algebra ([lq ) has a finite-dimensional represen- 
tation, of dimension 2j + 1, 

R 3 = { V \ V \ ...,^}. (16) 



The precise relation between the parameter j appearing in fll4]) and the QES pa- 
rameter J will be discussed in detail in the following. 
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It is convenient to collect in one place the inverse relations, connecting rf to 
the generators (|14]). These relations are 



and 



v 2 d 2 v 



_ T + T - 

T°T- + 



V 2 d v 
r)d v 



-(3j-l)T + + 2j(2j-l) V . 
+ 2jT° + 2j 2 , 



-T + + 2jt]. 
T°+J, 



(17) 



The matrix realization of the algebra (|15|) for semi- integer j, which will be useful 
for our purposes, is well-known: 












1 










2 




















2j 






' 





. 


. 


" 




2j 





. 


. 





T~ = 





2J-1 


. 


. 













. 


. 













. 


. 1 


_ 



j 
j-l 
j-2 






























-3 + 1 
-j 



(19) 



As was explained in Sect. 2, the algebraic sector of the periodic Lame system ([!]) 
consists of 2 J + 1 eigenvalues, which are the edges of the allowed bands or gaps. To 
relate the parameter j in (]14|) to the QES parameter J in (JJ) and (|3]) we need to 
distinguish between J even and J odd. 



3.1 J even 

For what follows it is convenient to introduce the operator 

H = + 4j(4j + 1) z/sn 2 (0|z/) - 2j(4j + 1) . (20) 

In this section we consider four si (2) algebraizations yielding the band boundaries 
for even values of J. If J is even, then J/4 is semi-integer. As we will see, in fact, in 
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this case the full algebraic problem is split into four completely disconnected distinct 
problems as follows: 

2J + l = {2(£) + l} + 3x{ 2 (^) + l}. (21) 

We will consider four distinct algebraizations, one with j = J/4 and three with 
j = 4 — ~. The first two groups in Eq. ( pT|) correspond to symmetric wave functions 
while the third and the fourth to antisymmetric wave functions. 



3.1.1 j = J/4, periodic wave function 

After the change of variables <p ~^ the Hamiltonian on the left-hand side of ([[]) 
takes the form 



H = H (j) = -4 (1 - u)tj - (1 - 2v)rf - vrf 



-2 



1 _ v ) _ 2(1 - 2u)rj - Ivrf d v + +l) u (l-r t )- 2j(4j + 1) 



= -4 v T+T° - 4(1 - 2u)T + T- - 4(1 - u)T°T~ 

-2 v (6j + 1)T+ + 4 (2j + 1)(1 - 2z/)T° - 2 (2j + 1)(1 - v)T~ 

+2j(l-2z/), (22) 

where j = j ■ The operator H(j) defined here is a basic element of the construction 
to be presented below. 

It is clear that the Hamiltonian H = H(j = ~) is purely algebraic; it acts on 
\&(z/) = P 2 j( z/ ), where P 2j - is a generic notation for a polynomial of degree 2j. There 
are 2j + 1 eigenfunctions of the type ^(v) = P2j(v). These eigenfunctions, ^(<p) = 
P2j(v ((/))), are periodic in <fi. Moreover, 2j + 1 eigenvalues are most conveniently 
found from the matrix realization (0), where j is taken as J/4; the dimension of 
the matrices T ±,Q is (2j + 1) x (2j + 1). This procedure gives the first contribution 
to the decomposition in (pip . 

3.1.2 j — periodic wave function 

It is known in the literature that the si (2) algebraization of the problem (|I]) admits 
more than one solution for the quasiphase e~ a . In Sect. [34M| the quasiphase was 
trivial, e~ a = 1. Now we choose another solution, e~ a = Jr)(l — 77). In other words, 



*( V ) = ^ V (l- V )i,( V ). (23) 
The quasigauge transformed Hamiltonian 

Hue, = e a He~ a (24) 



10 



acts on ip{rj). Then 

H«& = HQ) - 8uT + + 4(2 - 3z/)T° - 4(1 - v)T~ + {l + v + Av j } , (25) 
where 

- J 1 

J = 7 - 5 (26) 
and the generators T '° on the right-hand side of Eq. (p5"D, including those in H(j), 



are in the representation 2j + l (i.e., the matrices ( |l9l) of dimension (2j + l) x (2j+l)). 

We have to explain why the wave function ^(0) is periodic in this case. The 
situation is slightly more subtle than that in Sec. |3.1.1| . The eigenf unctions have 
the form 



* = ^ V (1- V )P 2 .( V ), (27) 

where 77 = r?(0). Inside each period the expression under the square root touches 
zero twice; at these points care should be taken of the branches of the square root. 
Needless to say that \l/(0) must be a smooth function of 0. Let us examine what 
happens, for instance, at near zero. It this point 1 — rj = 2 and \/l — rj must be 
understood as 0; the square root is positive at > and negative at < 0. This 
means that in the plane 77 we pass from one branch of the square root to another. 
This introduces a change of sign. Since this change happens twice inside the 
period, the wave function \l/(0) we deal with is periodic. If the change occurred 
once, the wave function would be antiperiodic. This is case for two algebraizations 
considered in the next subsection. 

3.1.3 j = '^-^1 antiperiodic wave functions 

There are two more solutions for the quasiphase e~ a , namely 

* = sj{l-rj)[{l-v) + vrj\ i>( V ), 

H«a» = #(j)-8z/T + + 4(l-2z/)T°-2(l-2z/)(l + 2]), (28) 

and 

* = ^r][(l-u) + u V } fa), 

H« G » = HQ) - Su T+ + 4 (1 - 3u) T° - 4 (1 - v) T~ 

- 2(l + 2j) + z/(l + 4j). (29) 

In both cases the expression under the square root touches zero once inside each 
period. At the point where it occurs one passes from one branch of the square root to 
another; correspondingly \l/(0) in the cases at hand is antiperiodic. The generators 
T '° in Eqs. (^Sj), ( P9"D are in the representation 2j + 1, while the eigenf unctions 
ip(r]) are polynomials of rj of degree 2j. This concludes our consideration of the band 
boundaries for even values of J. 
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3.2 J odd 

For what follows it is convenient to introduce the operator 
H(j) = -4 vT + T°- 4(1 -2v)T + T~ -4(1- u)T°T~ 

- 6u (2j + 1) T+ - 4(-l - 2j + 3u + Aju) T° - 2 (1 - u) (1 + 2 j) T" 

+ v-l-2j. (30) 

In this section we will consider four sl(2) algebraizations yielding the band edges 
for odd values of J. If J is odd, then (J — l)/4 is semi- integer. As we will see, in 
this case the full algebraic problem is split in four disconnected distinct problems of 
the following dimensions: 



4 J J I V 4 

The first two groups in Eq. ([H]) correspond to symmetric wave functions while the 
third and the fourth to antisymmetric ones. 

3.2.1 j = (J — l)/4, periodic wave function 

The wave functions ^f(4>) are periodic; they have the structure 



= yj[(l -U) + H ^(7?) = yJ[{\- V )+v n ] P 2J ( V ) . (32) 

The quasigauge-transformed Hamiltonian acting on ip(r]) has the form 

H"G" = H(j) , (33) 
with the generators acting in the representation of dimension 2j + 1. 

3.2.2 j = (J — 3)/4, periodic wave function 

There is another solution for the quasiphase leading to a periodic wave function, 



= Jri(l-ri) [(l-v) + vri] P 2 -( V ) , (34) 

where P 2 ](v) is a polynomial of degree 2j, 

~. _ . 1 _ J-3 
J - J - 2 - 4 

In this case 

ffoQ, = ft(j) -8 z/T + + 4(2-3 z/)T°- 4(1- z/)T~ 

+ 3z/ + 4z/j-l, (35) 
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3.2.3 j = (J — l)/4, antiperiodic wave functions 

The quasigauge factors emerging in two extra algebraizations leading to the an- 
tiperiodic are 



e a = ^fr] and y 1 — r\ . 
This implies the following wave functions and quasigauge-transformed Hamiltonians: 

* = Vv 4>(v) , 

H«g» = H(j) + (1-^)(1+4j) + 4(1-z/)T -4(1-z/)T-, (36) 

and 



H«Q" = H{j) + (l + 4j)+AT°, (37) 
where ip(r)) are polynomials of degree 2j. 

4 Particular examples 

It is instructive to consider a few simple examples. This consideration will help us 
to establish a general pattern of the band edge levels. 

4.1 J = 2, five band-edge levels 

The band (gap) edges in this case have been already presented in Eq. (^). According 
to the results of Sec. pO, the QES sector consists of one doublet (j = 1/2) and three 



singlets (j = 0). The doublet is obtained from Eq. fl22[) ; the corresponding energy 
eigenvalues are E\ and E 5 . Three singlets (E 2 , E s and E±) are obtained from Eqs. 
(p5f), (p8|) and (^9|). The doublet energies i?i and £5 are dual to each other; as far 
as the singlets are concerned E4 is dual to E2 (E4 corresponds to a periodic wave 
function while E2 to an antiperiodic wave function); and E3 is self-dual, see Fig. |2j. 

4.2 J = 3, seven band-edge levels 



According to the results of Sec. p72| the QES sector consists of three doublets 
(j = 1/2) and one singlet (j = 0). The doublet levels are obtained from Eqs. (|33l). 



and (|37|). 

£ 5 = -4 + 5z/ + 2Vl - z/ + 4z/ 2 



= -4 + 5^-2^1-^ + 4^, (38) 
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and 

E 6 = -1 + 2u + 2^4 - v + v 2 , 

E 2 = -1 + 2Z/-2V4-Z/ + Z/ 2 , (39) 

and 

E 7 = -1 + 5^ + 2^4-7^ + 4^, 

E 3 = -1 + 5z/- 2V4- 7z/ + 4z/ 2 . (40) 

The singlet is obtained from Eq. (|35"D, 

£ 4 = -2(1 - 2i/) . (41) 

The doublet energy eigenvalues E 2 and E% are dual to each other, E\ and are 
dual, and so are E5 and E 3 . The singlet energy eigenvalue E4 is self-dual, see Fig. 



ENERGY BANDS, J=3 




Figure 4: The energy bands for the Lame system (|l|), as a function of v. The 
plot given is for J = 3. The first (from the bottom) and the third allowed bands 
(unshaded) are of the p-ap type, while the second and the fourth ones of the ap-p 
type, see text. 



4.3 J = 4, nine band-edge levels 

According to the results of Sec. |3Tl] , the QES sector consists of one triplet (j = 1) 
and three doublets (j = 1/2). The triplet is obtained from Eq. ([22]), 

-^g 8\/l3 

#1,5,9 = 5- (1 - 2v) 5— Vl - V + V 2 COS [<J lj5i9 
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1 3^ V144 - 432z/ + 2089z/ 2 - 3458z/ 3 + 2089^ 4 - 432z/ 5 + 144z/ 6 

- arcsin = ■ 

3 26V^3(1 -p + u 2 f/ 2 



7T 

3 



7T 



5, = -- 



So 



7T . 



(42) 



At v < 1/2, the arcsine on the right-hand side is in the first quadrant, and is defined 
in a standard way. At v > 1/2 it must be defined as a smooth analytic continuation. 
The doublet levels are obtained from Eqs. fl25D , (|28|) and (^J), 



and 



and 



£ 4 



E 3 
E 7 



hv - 2V9 - 9u + 4u 2 



hv + 2V9 - 9u + 4u 2 , 



-5 + Wu - 2V4-9z/ + 9i 



-5 + lOu + 2V4 -9u + 9z/ 2 



(43) 



(44) 



E 2 

Er 



-b + hv- 2V4 + v + 4z/2 



-5 + 5z/ + 2^4 + z/ + 4z/ 2 . 



(45) 



The triplet levels ( f4~2"|) are dual, and so are the doublet levels E s and E 7 . Moreover, 
the two remaining doublets are dual to each other, namely £4 is dual to Eq, while 
E s to E 2 , see Fig. §. 



4.4 J = 5, eleven band-edge levels 



According to the results of Sec. |3.2| the QES sector consists three triplets (j = 1) 
and one doublet (j = 1/2). 

The triplet levels are obtained from Eqs. (|33"D , ( j36|) and (^). It is not difficult 
to obtain that 



£1,5,9 = o (~ 5 + 7zy ) - o V13 - 13z/ + 28z/2 cos £1,5,9 



1 9 v / 3\/16-48z/ + 461z/ 2 
- arcsin 



842z/ 3 + 541z/ 4 - 128i/ 5 + 256z/ 6 



2(13 - 13i/ + 28z/2) 3 /2 



7T 



7T 

3 ' 



7T . 



(46) 



This triplet is dual to that of Eq. (|4"8). As in Eq. (0), the arcsine on the right- 
hand side is defined in the standard way when it is in the first quadrant, (i.e. at 
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ENERGY BANDS, J=4 




Figure 5: The energy bands for the Lame system ([[]), as a function of v. This is 
the plot for J = 4. The first (the lowest), the third and the fifth allowed bands 
(unshaded) are of the p-ap type, while the second and the fourth of the ap-p type, 
see text. 

v < = 0.371...) while at larger v the arcsine is understood as a smooth analytic 
continuation 4 . 
Furthermore, 



10 8 
£2,6,10 = -— (1 - 1v) - - V28 - Yiv + 13^ 2 cos 
o o 



-'2,6,10 



1 9^^256 - 128i/ + 541z/ 2 - 842z/ 3 + 461z/ 4 - 48z/ 5 + 16z/ 6 

arcsm ; ——7- 

3 2 (28 - Yiv + 13z/ 2 ) 3 / 2 

7T 7T 

<5 2 = g, ^ 6 = - 3' *io = tt. (47) 

This triplet of levels is dual to itself. 
Finally, the third triplet is 



£3,7,11 = o (~ 2 + 7z/ ) ~ o V28 - 43z/ + 28z/ 2 cos 



3,7,11 



1 . 9^^256 - 1408z/ + 3741i/ 2 - 5162z/ 3 + 3741z/ 4 - 1408z/ 5 + 256z/ 6 

— clXCSlIl ; —7-777- 

3 2(28 - 43z/ + 28z/ 2 ) 3 / 2 



Eqs. @ and @ ^* = 1/2 and 0.629.., respectively. 
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5 3 =g, ^"3' 5 H =7r - (48) 

It is dual to the triplet (|46|) . 

The doublet levels are obtained from Eq. (f^), 

£ 4 = —5 + 10i/ — 6Vl — 1/ + z/ 2 , 

E 8 = -5 + 10Z/ + 6V1 - v + v 2 ■ (49) 

These two eigenvalues are dual to each other. The overall band structure for J = 5 
is presented in Fig. ||. 



ENERGY BANDS, J=5 




0.2 0.4 0.6 0.8 1 



Figure 6: The energy bands for the Lame system (^), as a function of v. This is 
the plot for J = 5. The first (the lowest), the third and the fifth allowed bands 
(unshaded) are of the p-ap type, while the second, the fourth and the sixth of the 
ap-p type, see text. 



4.5 The general pattern 

The remarkable features of duality and energy-reflection symmetry of the spectral 
problem ([]]) reveal themselves in a transparent manner in Figs. |2], f|, [| and || The 
allowed bands of the elliptic potential ^sn 2 (0|z/) present (up to a sign) the forbidden 
bands (gaps) of the potential (1 — v) sn 2 (0|l — z/), and vice versa. A nice illustration 
for duality of this type was given long ago by M. Escher. Figure |7] shows a fragment 
of Escher's "Sky and Water." Here the gaps between the birds are fishes, while the 
gaps between the fishes are birds. 
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Figure 7: An illustration of duality: a fragment of Escher's "Sky and Water." 



The last allowed band extends up to E = oo. The lower boundary of this band 
corresponds to periodic wave functions for even values of J, and to anti-periodic 
ones for odd values of J. The p-ap and ap-p allowed bands alternate, and so do the 
gaps. At v = 1/2, the energy bands are self-dual. This means that for each band 
(gap) edge with the energy E there is one with the energy — E. 

As v — > 1, the widths of the allowed bands (except the last one) tend to zero. 
In fact, these bands shrink to the discrete bound states of the Poschl- Teller system. 
The lowest band is at E = — J(J — 1)/2, the last band starts at E = J(J + l)/2. On 
the contrary, the gap widths decrease. The width of the lowest gap is AE = 2J — 1, 
and then it decreases: AE = 2J — 3, 2 J — 5, 1. 

As v — ► 0, the gap widths tend to zero, while the allowed band widths grow. The 
lowest allowed band is at -J(J+ l)/2 < E < -J(J + l)/2-l, its width is AE = 1, 
the next band is at -J(J + l)/2 - 1 < E < -J(J + l)/2 - 4, its width is AE = 3, 
and so on. The last (infinite-width) allowed band starts at E — J(J — l)/2. The 
bottom edge of the last (infinite-width) allowed band is of the p type for even J and 
ap type for odd J. 

In each (2j + l)-plet (j > 0) of the sl(2) representation (Sect. |3|) the energy level 
pattern is as follows: 

Ek , -Efco+4 , E ko+8 . 

The multiplet which includes Ej + i is special: it is always self-dual. The (J + l)-th 
energy eigenvalue E J+ i is dual to itself and passes through zero at v — 1/2. The 
dimension of this multiplet is 2 [J/4] + 1, where [...] denotes the entire part. 

The pattern described above is readily understood from the representation @ 
and/or the s/(2)-based representations considered in Sect. §. 
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5 Energy reflection symmetry at the self-dual 
point v = 1/2. 

When v — |, the Lame system is self-dual and so has an energy- reflection (ER) sym- 
metry. This symmetry implies that each eigenlevel of energy E is accompanied by a 
level of energy —E, the corresponding wave functions being related in a well-defined 
manner. A class of ER symmetric QES problems was constructed in [[/]. Although 
the construction of Ref. guarantees ER-symmetric spectra, by no means does 
it present a sufficient condition for the ER symmetry. The spectral problem (U), 
with v = 1/2, is in fact an expansion of the class of the ER-symmetric problems 
(the mechanism leading to this expansion, the possibility of different algebraizations 
of one and the same spectral problem due to the existence of distinct solutions for 
the quasi-gauge, was overlooked in 0). The focus of this section is the ER symme- 
try properties of periodic elliptic potentials from the standpoint of QES, and their 
connection with the Lame system ([[]). 

To begin, we observe that at v = 1/2 the function sn 2 (0|i/) is related to the 
Weierstrass function V with the invariants W3\ 



22 = 4, (?3 = 0; (50) 

namely, 

£^ {m ^ w _ l}s _£(^)_ 7 _L_ Y , (51) 

where 



•p ( _2 1 

' \V2 2 



_ 2^(5/4) 

r " r(3/4) 

is the period of the Weierstrass function V(x; gi = 4, g 3 = 0). The invariants #2,3 
will be suppressed below. This relation implies, in turn, that at v = 1/2 the spectral 
problem ([!]) can be written as 



1 d 2 J(J + 1) 1 



*(x) = £?*(x), (53) 



2 cfe 2 2 p _ 1) 

where the following change of variables is made: 

T^ x - (54) 

Note the occurrence of the factor 1/2 in the kinetic term. Weierstrass function 
related potentials first surfaced in the context of the Lie-algebraic construction in 
P], [TD| , and are well known in the general study of finite-gap potentials [T3| . 

In the general case of sn(^|z/) the parallelogram of periods is a rectangle; the 
v — 1/2 case is special: the parallelogram of periods is a square. This is called the 
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" lemniscate" case [13[]. The auxiliary variable ?y(0) 
becomes 

rj(x) 



1 



V(x) 



1 - sn 2 (0|z/) (see Eq. ([TJ)) 

(55) 



The analogue of Eq. ( ]I2"D becomes 



4(77 - r/ 3 



(56) 



A crucial relation 



??(ia;) = —r)(x) 



(57) 



follows from the properties of the Weierstrass function with the invariants (0). For 
real x the function rj(x) varies in the interval [0,1]; it is doubly periodic in the 
complex plane (along the real and imaginary axes), with periods r and ir. An 
example of the ER-symmetric potential (|53|) (J = 5), with the corresponding band 
structure, is shown in Fig. H. 



J=5, ER-SYMMETRIC 
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10 - 
5 - 




Figure 8: The ER-symmetric potential (|53|) versus x for J = 5. The gaps are shaded 
while the allowed bands are unshaded. Notice the ER symmetry which interchanges 
bands and gaps under E —> —E. 



5.1 J even, self-dual level multiplets 

Two solutions of this type for the band edges (i.e. individually ER-symmetric) were 
found some time ago in Ref. |7J (see Eq. (17) with v — and 1/2). The first one 
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(periodic eigenfunctions, cf. Eq. (^2|)) is 

3 = ~ v V(x) = P 2j (r ] (x)), (58) 
H = _ 2 (r j - V 3 )d 2 ri + (3 V 2 -l)d ri -2j(4j + l) V 

= -2T+T - 2T°T~ - (67 + 1)T + - (2j + 1)T~ . (59) 
The second solution (anti-periodic eigenfunctions, cf. Eq. Q2"8])) is 

j = = ^1 - [^)] 2 P$ , 

if« G » = -2( V - V 3 )d 2 TI + (7r } 2 -l)d ri -2j(4j + 5)r 1 

= -2T+T - 2T°T~ - (6j + 5)T+ - (2j + 1)T _ . (60) 

Note that the ER symmetry of each of the two level multiplets above follows from 
the fact that 

H -> -H 

under the replacement x —> ix, rj — > —77. The corresponding wave functions are 
obtained from one another by the replacement i] — > —rj. This is the mechanism 
discussed in |7j. In the example of Sec. [D^, J = 4, the band edges under discussion 
are 

{2^, 0, -2^13} , {V7, -V7}. 

5.2 J even, a level multiplet of the p type dual to that of 
the ap type 

Our assertion is: there exists another mode in which the ER symmetry can be 
realized in the problem (|I|). It does not fall into the category specified by Eqs. 
(6) and (8) in Ref. J7J. We consider two distinct quasigauge transformations, both 
leading to Lie-algebraic H«g" 's, such that the quasiphase factors e~ a are not invariant 
under 77 — > —T), but, rather, are transformed one into another. This will lead to a 
wider class of ER-symmetric problems than that considered in Ref. ]7) . The pair of 
conjugate Hamiltonians is 

j = 
+ 
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J-2 



y/ri(x) {1 - r](x)} P 2] (r](x)) 



-2T + T° - 2T°T~ 



2(3j + -)T + -(2j + 3)T- 



2T°+ [2j + 



(61) 



(periodic eigenfunctions, cf. Eq. (^5|)), and 

J-2 



j = —— , V(x) = ^rj{x) {l + r]{x)} P 2] (r](x)) 



H« G » = -2T+T - 2T°T- - 2 ( 3j + - ) T+ - (2j + 3) T 



- 2T° - (2j + |) , (62) 

(anti-periodic eigenfunctions, cf. Eq. (|29|)). 

The terms containing T° would be forbidden by the ansatz considered in Ref . [|7| . 
However, it is clear that the Hamiltonians (^1|) and ( |62| ) swap under the substitution 
i] — > —T]. In the example of Sec. £D| J = 4, the band edges under discussion are 



j- + V22, - - v^j , periodic ; j-- - y/22 , -- + v^j , anti-periodic 
5.3 J odd, self-dual level multiplets 

Two solutions of the self- conjugate type for the band edges are as follows: 

J- 1 



j = — — , V(x) = y/r^r) P 2j (r](x)) 



= -2T+T - 2T°T- - 2 (3j + T+ - (2j + 3) T~ , (63) 

(anti-periodic eigenfunctions, cf. Eq. (|36|) ), and 

J-3 



j = V(x) = yJ V (x) [1 - (v(x)) 2 } P 2] (r)(x)) , 

# lIG „ = -2T+T - 2T°T~ - 2 ^3j + T + - (2] + 3) T - . (64) 

(periodic eigenfunctions, cf. Eq. (^)). These two solutions fall into the category 
treated in Ref. [0. In the example of Sec. [14] , J = 5, the band edges under 
discussion are 



{2V33, 0, — 2V33}, anti-periodic; 
{3^, — 3v^3}, periodic. 

5.4 J odd, a level multiplet of the p type dual to that of the 
ap type 

The pair of conjugate Hamiltonians is 
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H«G» = -2T+T - 2T°T- - 2 (tlj + T + - (2j + 1)T~ 

- 2T° - (2j + ±) , (65) 

(periodic eigenfunctions, cf. Eq. (0)), and 

J- 1 



J = —4—, = V 1 - ^) P 2i (^)) , 

#« G » = -2T+T - 2T°T~ - 2 + -) T + - (2j + 1) T~ 

+ 2T° + + , (66) 



(anti-periodic eigenfunctions, cf. Eq. fl37|)). In the example of Sec. £Oj, J = 5, the 
band edges under discussion are 



-- - 4\/6cos U 159 - I + ^arcsinW^ ] « -11.49, -1.37, 5.36 }• , p; 



b - + iVdcos £3,7,11 " f + ^arcsinJ ^ ] = -5.36. 1 .37. 1 1 . 19 \ . a-]>. 



6 Weak coupling versus quasiclassical expansion 

In this section we describe how the duality transformation ([5]) connects various per- 
turbative and nonperturbative techniques, by relating information about states high 
up in the spectrum to states low down in the spectrum. We first consider approx- 
imate techniques for the locations of bands and gaps. Next we consider techniques 
for evaluating the widths of bands and gaps. The calculations of widths are sensitive 
to exponentially small contributions which are neglected in the calculations of the 
locations. 



6.1 Estimates of locations of bands and gaps 

Because of the duality transformation (|^), the location of a low-lying band in the 
spectrum is related to the location of a high-lying gap in the dual spectrum (i.e., the 
spectrum for the potential obtained by making the duality replacement v — > 1 — v). 
The location of the low-lying bands can be obtained by applying the weak coupling 
expansion. The location of a high-lying gaps can be obtained by applying the 
quasiclassical expansion. In both cases the quadratic Casimir J(J + 1) determines 
the expansion parameter. It is convenient to introduce a parameter k, 

K=yJj(J+l). (67) 
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Then 1/ k is the weak coupling constant of the perturbative expansion. Simultane- 
ously, 1/k plays the role of % in the quasiclassical expansion. 



6.1.1 Perturbation theory for location of lowest band 



In the limit J — > oo, the width of the lowest band becomes very narrow, so it makes 
sense to estimate the "location" of the band. In fact, as we will see in the next 
sections, the width shrinks exponentially fast, so we can estimate the location of the 
band to within exponential accuracy using elementary perturbation theory. That 
is, for large J, we can consider a single isolated well of the periodic Lame potential, 
and expand |TB[ near = 0, keeping quartic, sextic and higher order anharmonic 
terms, 



1 



sn 



+ 



2 + I3v + 2v 2 



45 



+ 



(68) 



Rescaling the coordinate, 



x 



r (ia/(J + 1))V4 
the Lame equation ([]]) becomes 

d 2 o v + 1 . 2 + \Zv + 2v 2 



+ X 



x 4 + 



x° + ... 



E + \k 2 t 



(69) 



dx 2 3k-s/u 45k 2 v k^jv 

Thus, the lowest energy level can be evaluated as a simple series in powers of 1/ac, 



En 



2^ u + 1 (l-Ais + v 2 ) 



2 k 2 



O 



1 



(70) 



6.1.2 Semiclassical method for location of highest gap 



As J — > oo, the location of the highest gap can be found by semiclassical techniques. 
First, note that for a given u, as J — > oo the highest gap lies above the top of the 
potential. Thus, the turning points lie off the real axis. For a periodic potential 



the gap edges occur when the discriminant [ IB] , |2T| takes values ±1. By WKB, the 
discriminant is 



A(£0=cos^ p h n S n (P) 



(71) 



where P is the period, and the S n (x) are the standard WKB functions (see e.g. 



S (x) = ^ y/Q(t)dt, 
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Si(x) 
S 2 (x) 

S s (x) 



--[logQ(x)]* , 



Q" 5(Q 



A2 



8Q 3 / 2 32 Q 5 / 2 



q" + 5(g') 2 



16Q 2 64 Q 3 



(72) 



and so on. Here 



Q(x) = E — V(x) . 

These functions S n (x) can be generated to any order by a simple recursion formula 
22]. The J-th gap occurs when the argument of the cosine in the discriminant ( |7T|) 



is Jti. This gives the location of the J-th gap, neglecting exponential corrections 
which correspond to the exponentially narrow width of the gap in the semiclassical 
J — ► oo limit. 

For the Lame system, we rescale the Lame equation ([!]) as 

" ^3?* + "-'<*">*= (f + 5) *• < 73 » 



and, thus, identify % in Eq. (fflj) as 



— <-> 7i . 



As a result, the condition for the occurrence of the J-th gap takes the form 



k £ —S n (2K) = Jn = nK(l- — 



n=0 



K 



1 



1 



1 



2 AC 8k 2 128 k 4 



(74) 



(75) 



where 2K is the period of the Lame potential, and where on the right-hand side we 
have expressed J in terms of the effective semiclassical expansion parameter 1/k. 
This relation ( |75| ) can be used to find an expansion for the energy of the J-th gap 
by expanding 



1 00 c 



(76) 



The expansion coefficients E\ are fixed by identifying terms on both sides of the 
expansion in ( fT5|) . For example, the leading terms clearly match because 



So (2K) 



2K E 1 



J™ y/l - z/sn 2 (0|z/) tfy + O (^j 



z/sn 2 (0|i/) d(p 



vr + 0(- 



(77) 
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The next-to-leading term on the left-hand-side of the expansion (F5|) comes from 
expanding S (2K): 



S (2K) = I™ Jl-usn 2 {(f>\u) + — d . . . 

JO V K 
7r£ l ^ ( 1 



. + Q - ■ (78) 

Matching the next-to-leading terms in (|75|), we find 

£l = -^\^. (79) 



This agrees with the next-to- leading term in the perturbative expansion (|70|), after 
making the duality transformations v — > 1 — v and E —>■ —E. 

Note that there is no contribution from Si, since Q(x) is periodic with period 
2K. (It is important here that the turning points are off the real axis, so that 
integrating along the real period we do not encounter any turning points, which 
would then require deformation of the integration contour in the complex plane 
PB| , P2"|). Indeed, for the same reason, none of the odd- indexed S n contributes to 
the discriminant. 

Similarly, successive orders in the energy expansion (^) follow by expanding 
(ff5|) to a given order in and matching to determine the expansion coefficients 
E£. In this way one finds 

1 o , 2-i/ (-2 + 2U+U 2 ) 

E = -k 2 -kVT^+ — + K - r J +■■■ ■ (80) 

2 4 16 Kyi — v 

Comparing with the perturbative expansion (f70|) we see that the semiclassical ex- 
pansion ( ^Ol) is indeed the dual of the perturbative expansion (ffOl), under the duality 
transformation v — > 1 — v and E — > —E. 



6.2 Estimates of widths of bands and gaps 

The semiclassical techniques used in the previous section were not sensitive to the 
exponentially small corrections needed to estimate the width of a low-lying band, or a 
high-lying gap. Yet, by the duality transformation the width of the lowest band 
is equal to the width of the highest gap, for the dual potential obtained by making 
the replacement u — > 1 — v. In this section we compute the widths of the lowest 
band and highest gap, using various approximate techniques, and compare with the 
exact duality result. This provides a link between perturbative and non-perturbative 
techniques that is much more sensitive than that discussed in the previous section 
for the locations of low-lying bands and high-lying gaps. 
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6.2.1 Algebraic approach for width of lowest band 



Since the band edge energies are given by the eigenvalues of the finite dimensional 
matrix H in (|]), the most direct way to evaluate the width of the lowest-lying band 
is to take the difference of the two smallest eigenvalues of H . For any given value of 
the elliptic parameter u, this involves finding the eigenvalues of a (2 J+ 1) x (2J + 1) 
matrix, which can be done with great precision. However, finding the eigenvalues as 
functions of v (see, for example, the plots in Figs. [| |], || and |6] for J = 2 and J = 5) 
rapidly becomes complicated as J increases. Nevertheless, from these expressions 
it is possible to deduce p(J the exact leading behavior, in the limit v — > 1, of the 



width of the lowest band, for any J 



A^gjraac = ^ ^ (1 - „)' ( 1 + -*,) + ...) (81) 



-.algebraic _ 8 JT(J + 1/2 ) M ,_\J f-, , J 

4^r(j) 

This clearly shows the exponentially narrow character of the lowest band in the 
v — > 1 limit. 



6.2.2 Tight-binding approximation for width of lowest band 



In the tight-binding approximation [ 2l| , one assumes the periodic wells are far sepa- 
rated, so that the periodic potential can be treated as a periodic sequence of "atomic" 
wells, each of which has a set of discrete bound levels. Small overlap effects broaden 
these discrete bound levels into bound bands. In the Lame case this approximation 
can be made very explicit due to the remarkable elliptic function identity: 



z/sn 



E' 
~K' 



( 71 



\2K' 



E 



sech 



/ 7T 



\2K' 



2nK)\ 



(82) 



Here E(u) is the complete elliptic integral of the second kind [|T^, D|, and 

E'M = E(l 



v) 



This identity ( |32"D shows that the periodic Lame potential can be written as a se- 
quence of periodically displaced (by the period 2K(u)) Poschl- Teller "atomic" wells 
(but the rescaling factor ^7 is non-obvious) . This "atomic" structure becomes clear 
graphically as v — > 1 (see Fig. |l|), but is in fact true for all v. Since each "atomic" 
well is a Poschl- Teller well, the normalized lowest energy bound state is 



W) 



y/¥r(j+i/2) 

\| 2K'T{J) 



sech 



/ 7T 



\2K' 



1 J 



In the tight-binding (TB) approximation, the energy of this lowest state broadens 
into a band of width fl2C 
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TB 

band 
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8jr(j + 1/2) 
4 J v^r(j) 



:i-u) 



j 



:i-v) + 



where in the second line we have kept dominant terms as v - 
that 

l + -(l-u) + ... 



1, and used the fact 
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This result (|84|) agrees precisely with the exact result ([81]) to this order in 1 — v. 
Therefore, for any J, the tight-binding approximation is good as v — > 1; i.e. as the 
separation between atomic wells becomes large. 



6.2.3 Instanton approximation for width of lowest band 



In the instanton approximation P% ], tunneling is suppressed because the barrier 
height is much greater than the ground state energy of any given isolated "atomic" 
well. For the Lame potential in ([!]), this means 



Remarkably, the instanton calculation for the Lame potential can be done in closed 
form P0| |, leading to 



A rpinstanton 
A ^band 



-2k 



i6) 



To compare this instanton result (|86|) with the algebraic result (]81|) we take z/ ■ 
and we take J to be large, in order to be in the semiclassical limit flS"5|). Then 



instanton 
band 



8J 3 / 2 
v^F4 J 



1 + 



J 



which agrees perfectly with the large J limit (using Stirling's formula) of the exact 
algebraic result (^T|). Thus, this example gives an analytic confirmation that the 
instanton calculation gives the correct leading large J behavior of the width of the 
lowest band, as v — > 1. For other values of v, the instanton formula dSBT) is also 



the correct leading large J result, but the comparison with A^!^ raic must be done 



numerically p(| for a given value of v since the large J asymptotic behavior of 
AE^^ ra,lc for arbitrary v is not analytically calculable. 



6.2.4 WKB approximation for width of lowest band 

We can also estimate the width of the lowest band using WKB. First, rescale the 
the Lame equation ([!]) as 



1 1 d ^ 1 2/ji \ t/j."\ ( 1 E 

2^^ + 2 sn ( « W) = U + 2^; 



= (88) 
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In this form we clearly see the small WKB parameter to be 

1 



'ft" 



so that the large J limit is indeed semiclassical. Then, the standard WKB [£5], |2C 
expression for the band width is 

AS WKB = ^ exp (- JL d( p fiiVM - S)) (89) 

where t.p. denotes the turning points. Thus, remembering to rescale the energy and 
using the first two terms in Eq. ([70]) for the energy of the lowest band, we get 



TT [ Jt.l 



'sn 2 (0|z/) \ . (90) 

Ky/V J 



Defining y = sn 2 (0|z/), we can now evaluate the integral in the semiclassical limit 
where "ft" = 1/(k^) < 1, 



1 V 1 ""*"/*' dy 1 l n f 1 + ^ 



-"ft" (-In "ft" + -ln(l-z/)-21n2-- 
\2 2 v ; 2 



0(("ft") 2 



(91) 



As a result, one finds [^(J from ( p0| 



AR 



WKB 
band 



instanton 
band 



(92) 



This factor of ye/ir has been known for a long time ^7j. (The instanton result 
for Ai^band is correct while Eq. fl9HD is off by this factor.) It has been found in 
many other comparisons between the instanton method and the WKB formula fl89"D 
n, m H II- It can be traced to a slightly crude matching of normalizations of 
wave functions [E7], pll p2|. A more precise WKB treatment leads, of course, to 



complete agreement between these two semiclassical methods, as has recently been 
emphasized explicitly for the Lame potential |33 . 



6.2.5 Algebraic approach for width of highest gap 

We now turn our attention to the width of the highest gap. Taking the difference of 
the two largest eigenvalues of the finite dimensional matrix H in (|3|), it is straight- 
forward to show that as v — > 0, for any J, this difference gives 



A palgebraic 
gap 
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(93) 
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which is the same as the algebraic expression (KTJ) for the width of the lowest band, 
with the duality replacement v — > 1 — v. 



This result illustrates a theorem due to Trubowitz which states that for a 
real analytic periodic potential, the gap widths shrink exponentially fast as one goes 
up in the spectrum. This in turn extends an earlier result of Hochstadt |35[|, which 



states that the gap widths go like O(^) if the potential is m times differentiable. 
For a real analytic potential, such as the Lame potential in ([[J), the gap widths 
shrink faster than any power of the gap label /, and, in fact, shrink exponentially. 
By duality, the low-lying bands are also exponentially narrow, as in QSll). 

6.2.6 Naive perturbative estimate for width of highest gap 

As v — > 0, for fixed J, the potential in ([!]) becomes weak. Thus, one should be able 
to solve this problem using perturbation theory. Near v = 0, the n-th gap can be 
considered as a splitting between the two degenerate free states = exp(±in</>) 
at E = n 2 — ~J(J + 1). A standard solid state physics approximation gives |2T[ the 
width of the n-th gap as (twice) the n-th Fourier component of the potential, 

AE^ p ^2\V n \. (94) 



The Jacobi elliptic function has Fourier decomposition [12 



EM 7T 2 00 n cos 

„_. M „ ) = 1 _|g + £ S _^. (95) 

Thus, for the highest (i.e., the J-th) gap, we would deduce 



K 2 sinh(^f) 



AE J » J(J+1) 



~ 8JV+1) [jq) , ^^0. (96) 

When J = 1, this agrees with the leading term in the exact algebraic result (|9~3|). 
But it does not agree when J > 2. This failure illustrates an important point - the 
formula (0) comes from first-order in perturbation theory. However, from (^) we 
see that the width of the highest gap is of J-th order in perturbation theory. So 
first-order perturbation theory is clearly not sufficient for J > 2. Indeed, to compare 
with the semi-classical (large J) results for the width of the lowest band, we see that 
we will have to be able to go to very high orders in perturbation theory. This is 
an interesting, and very direct, illustration of the well-known connection between 
non-perturbative physics and high orders of perturbation theory. 
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6.2.7 Perturbation theory to order J for width of highest gap 



It is generally very difficult to go to high orders in perturbation theory, even in 
quantum mechanics. For the Lame system (|]) we can exploit the algebraic relation 
to the finite-dimensional spectral problem ([3|). However, since H in (|3]) is a (2 J + 
1) x (2 J+l) matrix, the large J limit is still non-trivial. Since we are only interested 
in gap widths, in this section we neglect the constant term in H, and treat vJ 2 as 
a perturbation of the free matrix J|, which has eigenvalues 0, l 2 , 2 2 , . . . , J 2 . The 
non-zero eigenvalues of J 2 are two-fold degenerate, while the zero eigenvalue is 
nondegenerate. We are interested in the splitting of the two degenerate eigenvectors 
of J 2 with the highest eigenvalue, J 2 . Even though these two states are degenerate, 
they do not mix in perturbation theory, because of the form of the perturbation 
matrix J 2 . This is essentially because both and J 2 have a sub-block structure in 
which alternate rows and columns separate. For definiteness, we choose the following 
(2 J + 1) x (2 J + 1) representation: 
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diag ( J, J - 1, J - 2, . . . , - J + 2, - J + 1, - J) . 



(97) 



Then the matrix elements of J 2 between the orthonormal eigenvectors of J 2 have 
the simple sparse structure 
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where the nonzero entries in the top left corner are 
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(98) 



(99) 



The remaining diagonal entries appear in 2 x 2 diagonal blocks, diag(-D n , D n ), with 



D n = l -J(J+l)- l -{n + l) 2 



n 



1,2,..,J-1 



(100) 



Finally, the nonzero off-diagonal entries in (|98|) are 



O 



__ V / 2 J(J + 2)(J2 



n 



J -2. (101) 



In the basis used for (|9~8"D, the last two rows and columns refer to the degenerate 
states with the highest eigenvalue of Jj. We see that there is indeed no string of 
matrix elements connecting (J^) 2,7+1,27+1 with (Jy)2j,2j- Thus, these two degenerate 
states do not mix at any order of perturbation theory. 

Consider perturbing around the (2J)-th eigenstate of J 2 , with eigenvalue J 2 . 
The energy E 2 j of the perturbed system can be expressed implicitly as |36| 
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(102) 
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where E® is the n-th eigenvalue of the free system, and U m>n denotes the matrix 
element of the perturbation between the m-th and n-th free states. An analogous 
formula applies for E 2 j+i- 

^From the 2x2 sub-block structure of fl98|) , it is clear that for the terms in (|102 ) 
with strings of fewer than J matrix element factors, the shifts in E 2 j and E 2 j + \ 
are identical. Thus, these terms are irrelevant for the computation of the difference 
between E 2 j and E 2 j + i. So, we can concentrate solely on the J-th order term in 
( |102| ), which involves a product of exactly J matrix element factors. Also, since 
we are looking for the leading term at this order of perturbation theory, we can 
replace the true eigenvalue E 2 j or E 2 j + \ in the denominator by the free eigenvalue, 
which is J 2 . A further simplification follows from the fact that there is only one 
nonzero string of J matrix element factors beginning and ending with the index 2 J 
or 2 J + 1. Thus, to compute the difference E 2 j+i — E 2 j at J-th order, we simply 
have to take the difference between the string of matrix elements in each case, with 
the appropriate denominator factor. 

To proceed we need to specify whether J is even or odd. We consider J to be 
even (the case of J odd only requires simple modifications, and is left to the reader). 
With J even, there is in fact no nonzero string of J matrix elements beginning and 
ending at 2J + 1. However, there is such a string beginning and ending at 2 J, 



[^2J,2J-4^2J-4,2J-8 • • • ^8,4^4,l] 
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-j [J!(2J)(2J - 1)(2J— 2) ... (J + 2)( J + 1)] 



2(2J)! 
4 J 



(103) 



This, together with a factor of u J , gives the numerator of the J-th order term in 
the perturbation series ( |102j ). The denominator of the J-th order term, with E 2 j 
replaced by J 2 , is 



(J 2 - (J - 2) 2 )(J 2 - (J - 4) 2 ) ... (J 2 - 4) " J 2 = 4' 7 " 1 [(J - 1)!] 



(104) 



Therefore, the splitting between E 2 j + i and E 2 j is given, at order J in perturbation 
theory, by the ratio of (|T03| ) to ( |T0l|) , 
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(105) 



in complete agreement with the leading part of the exact algebraic result (|93|). 
This gives the leading v dependence of the width of the highest gap, for any J. 
This provides an explicit illustration of the connection between the non-perturbative 
results of the previous section, and high orders of perturbation theory. 
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6.2.8 WKB and "over the barrier tunneling" for width of highest gap 



As mentioned in section |6.1.2| , in the semiclassical limit, the highest gap lies above 



the top of the potential, so we need to use "over-the-barrier" WKB to determine the 
width of the highest gap. This means that the turning points lie off the real axis, 



and the width of the gap is given by the expression f25 | 

2 



A£™ = =- exp ( -21m J pdx) , (106) 

where p(<f>) = \J2(E — V(0)), and the integration is over a contour C beginning on 
the real <j) axis, then encircling one of the complex turning points, and returning to 
the real axis. The complex turning points lie on the vertical line through the point 
<p = K. Thus, it is convenient to change variables to 

(f) = K + %K' - i<f>' . (107) 

This is precisely the duality transformation change of variables (§) discussed in 
Sect. 0. Thus W% 01, 



E-V{(j)\u) = E + ^J(J+ 1) - J(J+ l)^sn 2 (0|z/) 

= E — -J {J + 1) + J {J + 1) (1 - v) sn 2 (0'|l - v) 

= V{<j)'\l-v)-{-E). (108) 

Therefore, since the location of the highest gap is approximately (see (Bjj)) 

1 



E = -J{J+l)-yj{l-u) J(J + 1), 

we see that the evaluation of the WKB integral in the exponent of (|106|) is identical 
to the WKB integral ( pT|) computed in Sect. |6.2.4| for the width of the lowest band, 



except for the replacement v — > 1 — v. Thus, the WKB estimate for the width of 
the highest gap is identical to the WKB estimate for the width of the lowest band, 
with the duality transformation v — > 1 — v. 



7 Conclusions 

In summary, we have shown that for the QES periodic Lame system ([!]), there is a 
duality symmetry that maps the spectrum into the energy-reflected spectrum of the 
dual Lame potential which has dual elliptic parameter 2/ = 1 — v. This means that 
bands or gaps high (low) in the spectrum are mapped into gaps or bands low (high) 
in the dual spectrum. The self-dual point, v — ~, of this duality transformation 
corresponds to an energy reflection symmetry of the self-dual potential. This also 
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provides an extension of the energy- reflection construction of M. Furthermore, this 
approach decomposes the (2 J + 1) x (2J + 1) algebraic spectral problem into four 
smaller algebraic problems (the precise details of this decomposition depend on 
whether J is odd or even). This decomposition represents a significant algebraic 
simplification in the large J limit. The large J limit is a semiclassical limit, and we 
have also shown in detail how the duality transformation relates the weak coupling 
(perturbative) and semiclassical (nonperturbative) sectors. Such a relation arises 
because the energy reflection aspect of the duality symmetry relates states high 
up in the spectrum to states low down in the spectrum. Interestingly, since the 
QES models discussed here are periodic, and hence have band spectra, this duality 
applies not just to the locations of the bands and gaps, but also to the widths of the 
bands and gaps. The calculations of such widths are sensitive to exponentially small 
contributions that are neglected in the previous calculations of locations of energy 
levels []5], [7|, . We were able to show that duality applies also to the calculation of 
these exponentially small effects. 

There are several directions in which it would be interesting to generalize this 
analysis. It would be interesting to find such a duality in models in more than one 
dimension. There is presumably a more geometrical interpretation of this duality, 
in terms of the associated genus J Riemann surface, and this may provide a useful 
alternative viewpoint which extends to other finite gap potentials. Finally, the 
Lame models possess other interesting transformation properties which remain to 
be explored from the QES perspective. For example, since 

J 2 X + vJl = uJ(J + 1)1 + (1 — u) (jl - j^jfj (109) 

the band edge energies for elliptic parameter v are related to those for elliptic param- 
eter — jz^i which reflects the behavior of the Jacobi function under such a parameter 
change. It would be interesting to investigate these more general modular transfor- 
mations in the context of quasi-exact solvability. 
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